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We use the next-to-next-to-leading order (NNLO) contributions to anomalous dimension governing the evolution 
of non-singlet quark distributions. We use the xF% data of the CCFR collaboration to obtain some unknown 
parameters which exist in the non-singlet quark distributions in the NNLO approximation. In the fitting 
procedure, Bernstein polynomial method is used. The results of valence quark distributions in the NNLO, are 
in good agreement with the available theoretical model. 



1. Introduction 



2. The theoretical background of the 
QCD analysis 
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The global partem analyses of deep inelastic scatter- 
ing (DIS) and the related hard scattering data are gen- 
erally performed at NLO order. Presently the next-to 
leading order (NLO) is the standard approximation 
for most important processes in QCD. 

The corresponding one- and two-loop split- 
ting functions have been known for a long time 
0, 1, I S, I S 0, I I M E3- The NNLO cor- 
rections should to be included in order to arrive at 
quantitatively reliable predictions for hard processes 
occurring at present and future high-energy colliders. 
These corrections are so far known only for the 
structure functions in the deep-inelastic scattering 
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Recently much effort has been invested in com- 
puting NNLO QCD corrections to a wide variety 
of partonic processes and therefore it is needed 
to generate parton distributions also at NNLO so 
that the theory can be applied in a consistent manner. 

S. Moch and et al. 0, E3 computed the higher or- 
der contributions up to three-loops splitting functions 
governing the evolution of unpolarized non-singlet 
quark densities in perturbative QCD. 

Du ring the recent years the interest to use CCFR 
data [18] for xF$ structure function in higher orders, 
based on orthogon al poly nomial expansion method 
has been increased E3 Hfl HE EH M 11 • 

In this paper we determine the flavor non- 
singlet parton distribution functions, xu v (x,Q 2 ) and 
xd v (x,Q 2 ), using the Bernstein polynomial approach 
up to NNLO level. This calculation is possible now, 
as the non-singlet anomalous dimensions in n-Moment 
space in three loops have been already introduced 

MM. 



In the NNLO approximation the deep inelastic co- 
efficient functions are known and also the anomalous 
dimensions in n-Moment space are available at this 
order p^. E3- Since in this paper we want to calcu- 
late the non-singlet parton distribution functions in 
the NNLO by using CCFR experimental data, we in- 
troduce the moments of non-singlet structure function 
up to three loops order. 

Let us define the Mellin moments for the i/N struc- 
ture function xF^,(x, Q 2 )'. 

M n {Q 2 )= [ x n - l Fz{x,Q 2 )dx . (1) 



The theoretical expression for these moments obey 
the following renormalization group equation 
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+7^(A- i )M n (Q 2 /^,A s (n 2 ))=0 

(2) 



where A s = a s / (47r) is the renormalization group cou- 
pling and is governed by the QCD /3-function as 
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P(A s ) = -2j2^s +2 



(3) 



The solution of Eq.© in the NNLO is given by 
1 



.4 



NNLO 



/^lnQnQVAf^) 
/9o 3 (lnQ 2 /A|^) 2 
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/3 2 ln(lng 2 /A|^)+/3 2 /3 -/3 2 ] . 



(4) 
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Notice that in the above the numerical expressions 
for /3 , [3\ and A are 

A) = 11 - 0.6667/ , 
fh = 102 - 12.6667/ , 

fh = 1428.50- 279.611/ + 6. 01852/ 2 , (5) 

where / denotes the number of active flavors. 

The solution of the renormalization group equa- 
tion for non-singlet structure function xF% can be pre- 
sented in the following form [221 ]: 



M n (Q 2 ) 
M n {Ql) 



exp 



A S (Q 2 ) An) 



MQl) 



1ns( x > 
P(x) 



dx 



C^ S (A S (Q 2 )) 



(6) 

where M n (Qo) is a phenomenological quantity related 
to the factorization scale dependent factor which we 
will parameterize in next section. is the anoma- 
lous function and has a perturbative expansion as 



NS 



(7) 



i>0 



At the NNLO the expression for the coefficient func- 
tion Cffg can be presented as [25| 



C%l(A s ) = 1 + C (1 \n)A s + C (2 \n)A 2 s 



(8) 



With the corresponding expansion of the anomalous 
dimensions, given by Eq. ||7J), the solution to the three 
loops evolution equation from Eq. ([6]), is as follows 



M 



NNLO 



(Q 2 ) 



MQ 2 ) 
MQl) 



l + [A s (Q 2 )-A s (Q 2 )} (3jL 
\A s (Q 2 )-A s {Qir Pl 



7<T 



m Kb. Po 



do 



+-[A 2 S {Q 2 )-A 2 S (Q 2 )\ 



A)" 72 Pi 11 ' Pi 

1 + C^A S (Q 2 ) + C^A 2 S (Q 2 )) V(n, Q 2 ) . (9) 
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Where V is the valence quark compositions as 

V{n, Ql) = u v {n, Q 2 ) + d v (n, Q 2 ) . (10) 

As we see in Mellin-n space the non-singlet (NS) parts 
of structure function in the NNLO approximation for 
example, i.e. A4^ NLO (Q 2 ), can be obtained from the 
corresponding Wilson coefficients and the non- 
singlet quark densities. 



In next section we will introduce the functional form 
of the valence quark distributions and we will param- 
eterize these distributions at the scale of Qq. 

By using the anomalous dimensions in one , two and 
three loops from [l|| and inserting them in Eq. 
the moment of non-singlet structure function in the 
NNLO as a function of n and Q 2 is available. 



3. Parametrization of the parton 
densities 

In this section we will discuss how we can determine 
the parton distribution at the input scale of Qq = 1 
GeV 2 . To start the parameterizations of the above 
mentioned parton distributions at the input scale of 
Qq we assume the following functional form 

xu v (x, Ql) = N u x a (l - x) b (l + cy/x + dx) , (11) 



xd v {x,Q 2 ) = ^(l-x) e xu v (x,Q 2 ) . (12) 

In the above the x a term controls the low- a; behavior 
parton densities, and (1 — x) h ' e large values of x. The 
remaining polynomial factor accounts for additional 
medium- a; values. Normalization constants N u and 
Nd are fixed by 



u v (x)dx 



d v (x)dx — 1 



(13) 



(14) 



The above normalizations are very effective to control 
unknown parameters in Eqs. (| 1 1 1 1 2[) via the fitting 
procedure. The five parameters with Aq^^ will 
be extracted by using the Bernstain polynomials 
approach. 

Using the valence quark distribution functions, the 
moments of u v (x, Qq) and d v (x, Qg) distributions can 
be easily calculated. Now by inserting the Mellin mo- 
ments of u v and d v valence quark in the Eq. (|10p . the 
function of V(n, Qq) involves some unknown parame- 
ters. 



4. Averaged structure functions 

Although it is relatively easy to compute the nth 
moment from the structure functions, the inverse pro- 
cess is not obvious. To do this inversion, we adopt a 
mathematically rigorous but easy method (2(| to in- 
vert the moments and retrieve the structure functions. 
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The method is based on the fact that for a given value 
of Q 2 , only a limited number of experimental points, 
covering a partial range of values of x are available. 
The method devised to deal with this situation is to 
take averages of the structure function weighted by 
suitable polynomials. These are the Bernstein poly- 
nomials which are defined by 



B n k{x) 



r(n + 2) 



Y(k + l)r(n — k + l) 



x k (l-x) n - k ■ n>k 



(15) 

Using the binomial expansion, the above equation can 
be written as 



B n , k(x) 



T(n + 2) 
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Now, we can compare theoretical predictions with the 
experimental results for the Bernstein averages, which 
are given by [H [1| 



F nk {Q 2 )= / dxB nk (x)F 3 (x,Q 2 ) 



(17) 



Therefore, the integral Eq. (|17p represents an aver- 
age of the function F^{x,Q 2 ) in the region x nk — 
^Ax nk <x<x n . k + TjAxnk- The key point is, the val- 
ues of F3 outside this interval have a small contribu- 
tion to the above integral, as B nk (x) tends to zero 
very quickly. In order to ensure the equivalence of 
the integral Eq. (fT7|) to the same integral in the range 
x\ = x nk - \Ax nk to x 2 = x nk + \Ax nk , we have 
to use the normalization factor, dxB n k(x) in the 
denominator of Eq. (|17|) which obviously is not equal 
to 1. So it can be written: 



F n k(Q 2 



f- 



k -lAx„L. 



dx B nk (x) F 3 (x,Q 2 ) 



-iAx n 



dx B nk (x) 



(18) 



By a suitable choice of n, k we manage to adjust to 
the region where the average is peaked around values 
which we have experimental data [l8l |. 

Substituting Eq. (TJJ|) in Eq. (JT7J) , it follows that the 
averages of F3 with B nk (x) as weight functions can be 
obtained in terms of odd and even moments, 



Fnk 



{n-k)\T{n + 2) 
T(k + l)r(n- fc+ 1) X 

n-k 

(19) 

We can only include a Bernstein average, F nk , if we 
have experimental points covering the whole range 
[xnk - \Ax nk ,x nk + ^Ax nk ] [HI]. This means that 



with the available experimental data we can only use 
the following 28 averages, including 



F, 



(cxp) 



(<2 2 )> ill (Q 2 )) -^4,2 (Q 2 ); 



f^I\q 2 ) 



Another restriction we assume here, is to ignore the 
effects of moments with high order n which do not 
strongly constrain the fits. To obtain these experimen- 
tal averages from CCFR data fisj ] . we fit xFs(x, Q 2 ) 
for each bin in Q 2 separately to the convenient phe- 
nomenological expression 



xFi 



phen) 



Ax B {l-x) c 



(20) 



This form ensures zero values for xF% at x = 0, and 
x = 1. In Table 1 we have presented the numerical 
values of A, B and C at Q 2 = 20, 31.6, 50.1, 79.4, 125.9 
GeV 2 . 



Q 2 (GeV 2 ) 


ABC 


20 


4.742 0.636 3.376 


31.6 


5.473 0.694 3.659 


50.1 


5.679 0.698 3.839 


79.4 


4.508 0.567 3.757 


125.9 


7.077 0.819 4.246 



Table 1 : Numerical values of fitting A,B,C 
parameters in Eq. (20) . 

Using Eq. (J^SJ with the fitted values of A,B and C, 
one can then compute F^ P \Q 2 ) in terms of Gamma 
functions. Some sample experimental Bernstein aver- 
ages are plotted in Fig.[TJin the higher approximations. 
The errors in the F^ P \Q 2 ) correspond to allowing 
the CCFR data for xF% to vary within the experimen- 
tal error bars, including the experimental systematic 
and statistical errors [18j • We have only included data 
for Q 2 >20GeV 2 , this has the merit of simplifying the 
analysis by avoiding evolution through flavor thresh- 
olds. 

Using Eq. (fT9)) . the 28 Bernstein averages F nk {Q 2 ) 
can be written in terms of odd and even moments. 
For instance: 



F 2 . 1 (Q 2 ) = 6 (M(2,Q 2 ) - M(3,Q 2 )) 



(21) 



The unknown parameters according to Eqs. () 1 II 1 2[) 

will be a,b,c,d,e and J^qcd- Thus, there are 6 pa- 
rameters for each order to be simultaneously fitted to 
the experimental F nk (Q 
subroutine MINUIT j| 
all the experimental data points and found an accept- 
able fit with minimum x 2 /dof = 74.772/134 = 0.558 
in the NNLO case with the standard error of order 
10~ 3 . The best fit is indicated by some sample 



QCD- 

x 

2 ) averages. Using the CERN 
|, we defined a global x 2 f° r 
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Figure 1: NNLO fit to Bernstein averages of xF$. 



curves in the Fig. [TJ The fitting parameters and the 
minimum % 2 values in each order are listed in Table 2. 



From Eqs. (|llll2[) . we are able now to determine 



the xu v and xd v at the scale of Qq in higher order 



corrections. In Fig. [2] we have plotted the NLO and 
NNLO approximation results of xu v and xd v at the 
input scale Qq = 1.0 GeV 2 (solid line) compared 
to the results obtained from NNLO analysis (left 
panels) and NLO analysis (right panels) by MRST 
(dashed-dotted line) [3l|] and A05 (dashed line) [321] . 
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Figure 2: The parton densities xu v and xd v at the input 
scale Qq — 1.0 GeV 2 (solid line) compared to results 
obtained from NNLO analysis by MRST (dashed line) 
[H[ and A05(dashed-dotted line) ^ . 





NNLO 


N u 
a 
b 
c 
d 


5.134 
0.830 
3.724 
0.040 
1.449 


N d 
e 


3.348 
1.460 


^qcd ' MeV 


230 


X 2 /ndf 


0.558 



Table 2: Parameter values of the NNLO non-singlet 
QCD fit at Qq = 1 GeV 2 . 

All of the non-singlet parton distribution functions 
in moment space for any order are now available, so 
we can use the inverse Mellin technics to obtain the 
Q 2 evolution of valance quark distributions which will 
be done in the next section. 



5. x dependent of valence quark 
densities 

In the previous section we parameterized the non- 
singlet parton distribution functions at input scale of 
Qq = 1 GeV 2 in the NNLO approximations by us- 
ing Bernstein averages method. To obtain the non- 
singlet parton distribution functions in x-space and 
for Q 2 > Qq GeV 2 we need to use the Q 2 -evolution 
in n-space. To obtain the x-dependence of parton dis- 
tributions from the n— dependent exact analytical so- 
lutions in the Mellin-moment space, one has to per- 
form a numerical integral in order to invert the Mellin- 
transformation [33j 

1 f°° 

f k (x,Q 2 ) = - dwx 
I" Jo 

Im[e i ' t 'x- c - wei ' k M k (n = c + we^,Q 2 )] , 

(22) 

where the contour of the integration lies on the right 
of all singularities of Mk{n = c + we l ^,Q 2 ) in the 
complex rt-plane. For all practical purposes one may 
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choose c ~ 1, 4> = 135° and an upper limit of integra- 
tion, for any Q 2 , of about 5 + 10/lna; -1 , instead of 
oo, which guarantees stable numerical results [H,[35|. 



6. Conclusion 

The QCD analysis is performed in NNLO based on 
Bernstein polynomial approach. We determine the 
valence quark densities in a wide range of x and Q 2 . 
Inserting the functions of q v (n, Q 2 ) for q = u, d in Eq. 
(22) we can obtain all valence distribution functions in 
fixed Q 2 and in z-space. In Fig. [3] we have presented 
the parton distribution xu v at some different values 
of Q 2 . These distributions were compared with some 
theoretical models (U HH . 

In Fig. [3]we have also presented the same distributions 
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Figure 3: The parton distribution xu v and xd v at some 
different values of Q 2 . The solid line is our model, 
dashed line is the MRST model [3l]], dashed-dotted line 
is the model from 1321. 



The QCD scale Aq^ 4 is determined together with 
the parameters of the parton distributions. In our fit 
results the value of Aq^ 4 and ol s {M%) at the NNLO 
analysis is 230 MeV and 0.1142 respectively. 
Complete details of this paper with calculation of LO, 
NLO and NNLO and also comparing them with to- 
gether is reported in Ref. [36]. 
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